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ABSTRACT 
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1 Introduction 



Integrable models in two- dimensions possess remarkable features. In particular the S- 
matrix for such a theory factorises into products of two-particle scattering matrices S^(9) 
(which gives the matrix element for the process a + b — > c + d) [28]. Under reasonable 
physical assumptions, it is possible to show that S must satisfy a number of rather simple 
algebraic equations. Solutions have been found to these equations for all affine Toda 
theories and these are postulated to give a non-perturbative expression for the scattering 
matrix of the theory [5, 10, 7]. 

These ideas have been extended to include integrable theories on a half-line [6, 13, 16]. 
In this case it is necessary to introduce another matrix K which describes the reflection 
of particles off the boundary. Again it is generally believed that all scattering amplitudes 
factorise into products of S and K. In addition the algebraic equations which were solved 
by S alone previously can be modified to include K. Recently a number of solutions to 
these equations for S and K have been found in the context of affine Toda theories [8, 9, 13, 
14, 19, 24]. However, many of these solutions make no reference to the boundary conditions 
which presumably need to be imposed to make sense of the theory. (The exceptions to this 
are [19] whose perturbative analysis is tied to Neumann boundary conditions and [8,9]). 
In fact, integrability places severe constraints on the boundary conditions one can impose; 
for simply-laced affine Toda theories there are only a finite number of possibilities [8, 9, 3]. 

Fortunately, there seems to be a straightforward way of analysing which solutions for 
S and K correspond to which boundary conditions . Whilst S is known to tend to unity 
in the classical limit, it seems that in general the reflection factor K does not. Thus if 
we know the classical reflection factors for a particular boundary condition, it should be 
possible to identify the corresponding quantum reflection factor by considering its classical 
limit [8]. 

The aim of this paper is to make some progress towards calculating classical scattering 
for Toda theories on a half-line. This problem naturally divides into two parts. First it is 
necessary to find the 'vacuum configuration' or lowest energy solution (which is presum- 
ably static) which satisfies a particular boundary condition. Then one solves the linear 
equations for infinitesimal perturbations about this solution still ensuring that the bound- 
ary conditions are satisfied. Far from the wall the solution consists of a superposition of 
incoming and outgoing waves. The relative phase between the waves is interpreted as the 
classical limit of the reflection factor K. Some examples of calculations along these lines 
are to be found in [8] where particular cases within cfip and were considered. 

The paper will be divided into five sections. In the next section we review affine 
Toda theories, and in particular the boundary conditions which can be imposed which are 
consistent with the (classical) integrability of the theory [3,4]. It is then shown that the 
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boundary conditions combine in a particularly neat way with the equations of motion to 
give linear equations that the tau functions must satisfy at the boundary. This is the key 
result of the paper which enables us to solve for the classical scattering solutions. 

In the third section we consider static solutions to the boundary conditions for the 
special case of Toda theories. For these theories a large number of tau functions can 
be constructed explicitly by analytically continuing the multi-soliton solutions of the imag- 
inary coupling theory [18]. Whilst it can be shown that these solutions are singular on 
the whole line, it is possible that all the singularities can be placed 'behind' the boundary 
for a theory on the half-line. By substituting this family of tau-functions into the equa- 
tion derived in the previous section we obtain a large number of possible solutions. The 
correspondence with particular boundary conditions is discussed. 

In section four, whilst remaining within the context of affl Toda theories, the analysis 
is extended to cope with the scattering perturbations around vacuum solutions. In this 
way we are able to derive an expression for the classical reflection matrix. It is shown that 
this satisfies the classical reflection bootstrap equation. 

We conclude with comments on our results and directions for future research. 

2 Classical affine Toda theory on a half-line 

In this section we briefly review some of the features of affine Toda theory on a half line 
that we shall need later on to establish notation. 

To each affine Kac- Moody algebra g we can associate an affine Toda theory [22]. For 
simplicity in this paper we shall restrict ourselves to simply-laced algebras, and for the 
most part to affl which are in some senses the simplest cases of all. The equations of 
motion for the affine Toda theory associated to g are given by 

V 

Tfl - 

d^<f>+-'£n i a i e (3aK < t > = 0. (2.1) 

Here we have used the notation that a 1 for % — 1 to r are the simple roots of g, the finite 
Lie algebra associated to g, and a = —ip where ip is the highest root of g (for untwisted 
algebras). Also we have defined r to be the rank of g and introduced m the mass parameter, 
and (3 the coupling constant of the theory. In this paper we shall consider theories for which 
(3 is real, and henceforth we set (5 = m = 1. The constants rii are given by the equation 

r 

a + E n ^ = ° ( 2 - 2 ) 

i=i 

and are sometimes referred to as 'marks'. 
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When the classical theory is considered on the half line — oo < x < 0, we need to 
supplement the equations of motion by boundary conditions. In this paper we shall restrict 
ourselves to conditions of the form 

d x (j> = F{(j>) (2.3) 

although more general possibilities can be considered [27,4]. An arbitrary choice of F((f)) 
will generically break the integrability of the system. It turns out that for simply-laced 
theories (with the exception of sinh-Gordon or affine Toda theory [21,26]) there are 
only a finite number of choices for F which preserve integrability [8,9,3]. These can be 
summarised as follows: 

(1) Either 

OA = 0; (2.4) 

these are 'free' or Neumann boundary conditions 

(2) Or 

d x( p = iTA l a^ t e aK ^ 2 (2.5) 

i=0 

where the A { = ±1. 

The first of these two possibilities is fairly easy to analyse completely at a classical 
level. The classical energy functional on the half-line is 

E = /^^Q^ + ^^^^n^e^-l)). (2.6) 

This is non-negative and the lowest energy configuration is clearly = which also satisfies 
the Neumann boundary conditions. If we now take the field <f>(x, t) = e(x, t) to be some in- 
finitesimal perturbation to this vacuum configuration we see that in a linear approximation 
we have 

d^e + Me = (2.7) 

where M = J2l=o n i° L% ® a% i s the mass matrix. The eigenvectors p a , a = 1 to r, of M are in 
one-to-one correspondence with the fundamental representations of the Lie algebra g and 
are interpreted as the basic particle-like excitations of the theory. The eigenvalue of p a is 
given by A a = m 2 a where m a is the mass of the corresponding particle. It is a remarkable 
fact that the set of masses form the lowest eigenvalue eigenvector of the Cartan matrix of 
g [5, 11, 12]. In terms of p the basic solutions of 2.7 are given by 

<f>{x, t) = e(x, t) = p a e E \^ x + K a e-^) (2.8) 

where E 2 —p 2 = m\. This is a superposition of left- and right-moving waves. The classical 
reflection factor is given by the phase factor K a which can be determined by substituting 
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the solution (2.8) into the boundary conditions (2.4). This yields the solution K a = 1. 
Classically, at any rate, this boundary condition seems rather uninteresting. There has 
been considerable progress in understanding the quantum case in [19], where the semi- 
classical reflection factor has been calculated using perturbative techniques. 

In this paper we shall be interested in boundary conditions of the form (2.5). It turns 
out to be particularly convenient to analyse this case in the language of tau-functions [18]. 
These are introduced as a particular parametrisation of the function 0; 

= -][> l Zn(r J ) (2.9) 

i=0 

Note that we have introduced r + 1 functions Tj to describe the r-component field and 
this is reflected in the freedom to scale all the tau-functions Tj — > f(x, t)rj without affecting 
the value of 0. This unphysical degree of freedom is fixed by demanding that r« satisfies 
the following particular form of the equations of motion 

nn - n 2 - rfr, + [r[f = (jjr^ - r 2 j n h (2.10) 

where is the incidence matrix of g given by 

/, ; ■2d; J ot-a? (2.11) 

Essentially, iy takes the value one if the nodes corresponding to i and j are connected on 
the extended Dynkin diagram, and vanishes otherwise. Substituting the form (2.9) into 
the boundary conditions (2.5), and taking the inner product with the fundamental weight 
Aj, we discover that the boundary conditions can be written in terms of the tau-functions 

as 

^ + v/^V"^ 2 = n t C (2.12) 



Ti 



where 



From this we see that 



C = (^ + A e a °^ /2 ). (2.13) 

T~0 



Hi 
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11^ =r^ ni C-^r (2-14) 



Substituting this into the equations of motion (2.10) we find 



in) 



2 



r'l + 2 ni CT[ - (ntC 2 - n^n = 0. (2.15) 
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This equation is the main result of this section. Immediately we see that it has two 
attractive features. Firstly the equations for T; L have decoupled so that each equation only 
involves 7$ for some %. Secondly, if we assume that the solution is static then the first two 
terms vanish and we are left with the linear equation 

t'I - 2mCrl + (n\C 2 - n t )n = (2.16) 

At this point we perhaps should remind the reader that this equation is deduced from the 
equations of motion and the boundary conditions, and that the latter are only valid at the 
boundary x — 0. Thus the above equation cannot be solved as a differential equation in x 
since it is only valid at the boundary, but it will prove a valuable tool in determining the 
subset of solutions of the equations of motion which satisfy one of the boundary conditions 
(2.5). One of the drawbacks of the equation (2.16) is that we have essentially squared the 
boundary conditions, so that we have lost the information contained in the A^. This we 
will have to recover by explicitly examining the proposed solutions at the boundary 

Perhaps even more surprisingly we can use a similar equation in the case of time- 
dependent perturbations to a static vacuum solution, i.e. in the situation where 

Ti(x, t) = (Ti(x)) vac + €i(x, t) (2.17) 

where €i(x,t) is some infinitesimal perturbation to the vacuum solution {ji) vac . In this 
case, the important point is that fj = 0(e), so discarding terms of 0(e 2 ) in (2.15) we again 
arrive at the linear equation 

f t - t'I + 2mCTi - {njC 2 - m)Ti = (2.18) 

We shall use this equation in section four to solve for classical scattering about the vacuum. 

One may be concerned that one has introduced a spurious constant C in equation 
(2.16). Remarkably C has a physical interpretation as proportional to the energy of the 
solution corresponding to Ti(x, t) on the half-line. The energy on the half-line for boundary 
conditions of type (2.5) is given by 

E = E bulk -j^lA^e^-^X^ (2.19) 

i=0 

where Ebuik is the bulk energy given by (2.6) and the second term is the contribution to 
the energy from the boundary compatible with the boundary conditions (2.5). It has been 
shown in [23] that for soliton solutions the energy density is a total derivative, so that the 
bulk energy can be written as a boundary contribution which in terms of tau-functions is 



E, 



bulk 



-2£ 

i=0 



(2.20) 



J x=— oo 
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In the sections that follow we shall be using tau-functions that tend to a constant as 
x — > — oo so the contribution to (2.20) come only from x — 0. Adding the contribution 
from the bulk and the boundary we find that 

E = - 2 j2(^- + A t ^r t e a ^/ 2 ) (2.21) 

r 

= -2Y,niC = -2hC (2.22) 

i=0 

where h is the Coxeter number associated to g. Thus, although initially it may have seemed 
that the appearance of C in the equation (2.16) was a drawback, it turns out that it is an 
added bonus, giving the energy of the solution that we are considering. This is important 
in trying to determine the 'vacuum' for a given boundary condition, since the vacuum 
is defined to be the static solution with the lowest energy which satisfies the boundary 
condition. 

3 Static vacuum configurations for affl Toda theories from 

analytically continued solitons 

In this section we shall confine ourselves to affine Toda theories based on the algebra affl . 
We have seen in the previous section that the requirement that a solution satisfies one of the 
integrable boundary conditions of the form (2.5) can be neatly expressed in terms of tau- 
functions. For the we have a particularly rich source of such tau-functions which can 
be constructed by analytically continuing multi-soliton solutions of the imaginary coupling 
theory [18]. 

At this point let us briefly review the nature of these solutions in the imaginary coupling 
theory. Let us reintroduce the coupling constant, so that we write 

P<P=-j2^Mri) (3.1) 

i=0 

The tau-functions for an iV-soliton solutions can be compactly written as 

1 1 / n \ 

Tj (x, t) = ]T . . . ]T exp ]T ^% + Yl W q lnA^ . (3.2) 
iui=0 hn=o \p=l i<p<q<N J 

Here we have introduced the notation $ p = a p (x — v p t) + £ p where v p is the velocity of 
the p-th soliton, £ p is a complex parameter whose real part and imaginary part are related 
respectively to the position and the topological charge of the p-th soliton. Also o v and v p 
are related by the mass-shell condition 

ol(l-v})=ml p (3.3) 
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where a p labels the species of soliton and m ap = 2 sin(^i). We define uo = exp(2ni / (n+1)) . 
The variables a p and v p are often conveniently parametrised in terms of the two-dimensional 
rapidity 6 P by putting 

Op = m ap cosh(6» p ) (3.4) 
a p v p = m ap sinh(6> p ) (3.5) 

The constants A^ apaq " > describe the interaction between the p-th. and q-ih solitons and are 
given as 

, , (a v — a a ) 2 — (a v v v — a a v a ) 2 — 4 sin 2 -?-r( cir, — a a ) 

j^(a p a q )^Q^ _ V P 1> \_P_P 9 9/ n+1 V V 1> ^ 

(o-p + a q ) 2 - (a p v p + (T g v q ) 2 - 4 sin 2 ^ ( Gp + a<? ) 



blI H 2 i+ 2(n+l) ^ blI H 2 i 2(n+l) j 
Sm V 2 j + 2(n+l) ^ Sm V2i 2(n+l) ^ 



(3.7) 



where 6 = ^ p — In this section we shall only be interested in static solitons so we shall 
take 6 P = or equivalent ly v p — 0, <7p = m ap . 

For single solitons the expression (3.2) reduces to 

tj = 1 + u ]a e^ (3.8) 

The topological charge of such a soliton of species a which is defined by 

Jim (0(x, t) - (f)(-x, t)) (3.9) 

can be shown to lie in the fundamental representation with highest weight A a where X a -a b = 
8 b a . It is therefore natural to associate each species of soliton with nodes on the Dynkin 
diagram of a n . By analogy with the representation theory of a n we refer to solitons of type 
a and type a = n+ 1 — a as conjugate. We shall (ab)use the notation p to denote a soliton 
whose type is congugate to that of soliton p. 

Now let us specialise to the case of interest, namely static solitons. Examining the 
expression (3.2) carefully we see that solitons obey a kind of Pauli-exclusion principle. We 
can only construct multi-soliton solutions whose constituent solitons have either different 
velocity and/or different species. If we attempt to consider two solitons of the same species 
and velocity, the interaction constant A^ aa \0) vanishes and we simply end up with one 
constituent soliton of that velocity and species at some different position. This places 
severe constraints on the possibilities for static soliton configurations; since the velocities 
of each constituent is identically zero, it follows that each constituent soliton must be of 
different species. Thus we can consider at most an n-soliton solution where each constituent 
soliton corresponds to a different node on the Dynkin diagram of a n . Indeed we believe 



7 



that this is the most general static solution whose energy-density tends to zero at spatial 
infinity The energy of a stationary soliton of type a on the whole line is simply given by 

E = ~(n + l)m a (3.10) 

Note that for imaginary (3 this is positive, as we might expect. The energy of N stationary 
solitons is simply given as the sum of energies of the constituent solitons E = J2 p E p . 

The above discussion has all been in the context of imaginary coupling Toda theory, 
where it is accepted that the field is allowed to be complex. We are interested in real 
coupling Toda theory where the field is is taken to be real. A little thought shows that 
for tau-functions of the type given in (3.2), reality of implies the reality of Tj. This can 
be ensured by insisting that each constituent soliton p is paired with a congugate soliton p 
with the position/topological charge variables related by £ p = (Cp)*- An exception to this 
rule occurs for n odd, where the soliton corresponding to the middle node of the dynkin 
diagram, i.e. a p = (n + l)/2, is unpaired (since it is its own conjugate), so we must 
take£ p to be real in this case. Thus the reality of the tau-functions closely corresponds 
to the representation theory of a n , where only the middle node corresponds to a real 
representation and the other fundamental representations must be taken in conjugate pairs 
if we wish to restrict ourselves to real representations. Similar remarks apply equally well 
to tau-functions associated with other simply-laced algebras 

Finally, let us note that the energy E for such solutions on the whole line given in (3.10) 
is negative for real [3. This may seem surprising in view of the manifestly positive energy 
density of the theory (2.6). This reflects the fact that all such solutions become singular 
somewhere on the real line, and this is perhaps why they are not generally considered in 
the real coupling theories. However, in the present context the possibility exists that all the 
singularities in lie in the 'unphysical' region x > 0, so that they are perfectly acceptable 
physically. 

3.1 TWO-SOLITON SOLUTIONS 

As a pedagogical introduction to a more general solution, we begin by considering the 
possible static two-soliton solutions that satisfy the boundary conditions (2.5). This solu- 
tion and some of its features have already been discussed in [15]. A two-soliton solution 
consisting of a soliton of type a and its anti-soliton of type a have tau-functions of the form 

n = 1 + 2dcos(x + ^-) e maX + A {aa) d 2 e 2maX (3.11) 
v n+1 v ' 
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Here we have set x = I m {0 an d d = exp(i?e(£)) in the previous notation. Also we 
calculate from setting to zero in (3.7) as 

A (™) = cos 2 (— ) = 1 - (3.12) 

In general multi-soliton tau-functions (3.2) can be split up into a sum of 'charge' sectors 
by writing 

Tj = Y, T ^ kj - (3-13) 

k 

The linearity of equation (2.16) and the fact that for a n , we have that n £ = 1 for all % 
implies that each of the T k separately satisfy the equation; 

(T k )" - 2C(T k )' + (C 2 -l)T k = at x = (3.14) 

In the two-soliton case there are three charge sectors: k — 0, k — ±a. Obviously T a = TL a , 
so there are essentially two independent equations; 

(C 2 - 1) + ((C - 2m a f - l)A (ad) d 2 = (3.15) 
((C-m a ) 2 - l)d = (3.16) 

The second of these equations implies that for a non-trivial solution 

C± = m a ±l, (3.17) 

and, substituting this into the first we find that 

d = - — - — -. (3.18) 

The energy of the two solutions are given by E± = —2hC± = — 2(n + l)(m p ± 1). 
Actually, it is clear that solutions should come in pairs. The reason is that if 4>(x, t) is 
a solution to the equations of motion and the boundary conditions, then so is the parity 
reversed solution <f)(—x,t), since it satisfies boundary conditions of the form (2.5) but 
with Ai — > —Ai. Note that if we sum the energies of the two solutions, the boundary 
contributions of each solution cancel, and we are left with a bulk contribution of <p(x, t) 
and (f>(—x, t) from — oo < x < which can be rewritten as simply the bulk energy of <p(x, t) 
on the whole line, and indeed we find that 

E + + £_ = -2(n+ l)m a (3.19) 

which is the energy of two static solitons of type a. 
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Before we pick the lower energy solution as the vacuum we should be careful to consider 
whether the singularities of the corresponding solution for <fi lie in the physical region x < 
or the unphysical region x > 0. Since we have established that the two solutions are parity 
conjugate, and we know that the solution is singular at some point on the whole line, it 
follows that at least one solution must be singular somewhere in the physical region. (It is 
possible that the singularity lies at x = which is energetically allowed. We discuss this 
case later). It turns out that the 'good' solution whose singularities lie in the region x > 
corresponds to the higher energy solution C_. For this solution the tau-function can be 
written 

(2 + m a ) + 4 cos(x + gf)e m °* + (2 - m a )e 2m "* 

2 + m a {6 - 2()) 

Singularities in occur when Tj vanishes for some j. Solving the equation Tj = yields 



emaX = -cos^+^^y/cos^ + ^-cos^ 



For a real solution for x this implies that 



2nja 2 na 



cos>(x + ^fj) > co^(— ) (3.22) 

cos( X + ^-) < 0. (3.23) 
Combining these two conditions gives 

_1 < C os(x + — ^) < - cos(^ T ) (3.24) 

It is easy to check that for any choice of x there will be some j for which this condition is 
satisfied, so that there will always be a singularity at some real value of x for some j. To 
see where these singularities lie consider the right hand side of (3.21) where the square-root 
is taken with the minus sign, since this is the smaller of the two solutions, and hence the 
more likely to yield negative values of x. A simple calculation shows that as a function of 
cos(x + ) the right hand side is a monotonically increasing function, so its lowest value 
is attained when cos(x + ~rf) — — 1- At this point the right hand side is equal to one, 
corresponding to a singularity at x — 0. For other values of cos(x + f^pf), the right hand 
side will be greater than one, and the singularities will be in the unphysical region x > 0. 

Having found a class of two soliton solutions which correspond to non-singular fields 
0, it only remains to ascertain precisely which boundary conditions each solution satisfies. 
Substituting x = into the expression (3.20), we find that the tau-function at this point 
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Figure 1: Boundary conditions for N — 1, n — 8 and a = 2 



can be written as 



cos 2 (f + ^) 



Thus we can write 



s . cos 2 (f + ^±i^)cos 2 (f + ^^) 
" C os 4 '*^^ " " ~ ' 



2 ~ n+l> 



In the boundary conditions (2.5), we have terms involving exp(ai-? ■ 0/2) which should be 
interpreted as the positive square root of the above quantity. The values of Ai are found 
to be given by 

A, = -sign (cos(- + — rr ~) cos(- + -^3-) J (3.27) 

In Figure 1. we have plotted cos(| + ^f), tabulated its sign, and tabulated the sign 
of the right hand side of (3.27). From this we can see that the Aj is generically negative 
except at near a point where cos(| + ^-) has a zero, where a pair of positive Aj are 
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produced. The number of positive pairs is related to the number of such zeroes which are 
in turn related to the value of a. One also observes that the expression (3.27) is invariant 
under 

x ~* x ~TTTi < 3 - 28 > 

j -. j + l (3.29) 

This implements the Z n symmetry of Toda theories, and cyclically permutes the 
boundary conditions Aj — > mo d n +i ■ 

3.2 MULTI-SOLITON SOLUTIONS 

In this section we shall generalise the two-soliton solution given above to multi-soliton 
solutions with arbitrarily many constituent solitons. Whilst this case is much more com- 
plicated, and we cannot guarantee finding a complete set of solutions to the basic equation 
(2.16), we shall give a wide class of explicit solutions. 

To begin with we shall restrict our attention to multi-soliton solutions with an even 
number of constituent solitons; that is, we shall assume that the soliton corresponding 
to the middle node of the Dynkin diagram for n odd does not feature. Thus our ansatz 
consists of N pairs of conjugate solitons which we shall label p and p respectively. The real 
tau-function can be written in the form 

N <2 

1 + 2 ]T d p cos( Xp + —?-)e m ^ + ... (3.30) 
p=i n + 1 

where the dots compactly represent the (many) interaction terms. As in the two-soliton 
case we can decompose the tau-function into sectors as in (3.13). As a further simplifica- 
tion we shall assume that n is large with respect to the charges of the solitons, or more 
particularly that the largest charge Qmax = Y^pO-p < (n + l)/2. With this restriction, the 
charges in (3.13) take values between Q max and — Q ma x, and only the term 

n(rf p e^ + w))( J] A^))exp(f> ap x) (3.31) 

p=l ^ ' \l<p<q<N J \p=0 / 

contributes to the charge Qmax- Substituting this into the equation (3.14), we immediately 
find 

N 

C± = ]Tm ap ±l (3.32) 

p=0 

and the energy is as usual — 2{n + 1)C±. Once more this pair of solutions are the parity 
inverses of one another. To determine the 'positions' of the constituent solitons we consider 
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the terms with charge Q p = Q m ax — CLp- Again we make the assumption the charges are 
chosen in such a way that only two terms in the Tj contribute to Tq p : namely the term 
containing all the solitons except the p-th, and the term containing all the solitons and the 
p-th conjugate soliton. More explicitly we have 



tq p = n (D r j<* + ^) n Aiaras) 



exp I ^2m ar x 



r^p \r<s:r,Sy£p J yr^p 



x ( 1 + d 2 p A pp J] (A {rp) A (rp) ) e 2m ^ x j (3.33) 



Substituting this into (3.14) we discover that 



{d p ) 2 A pp J] (A irp) A irp) ) = 2±ma " (3.34) 
r+p ^ =F m ap 



One can prove that 



so that 



A (rp) A (rp) = K m a r m a J 

(m ar +m ap Y 

, ttK+^)J_ {3M) 

r^p \ m a r - m a p \ ^ T m ap 

Here we have made a choice for the sign of d p , since we can absorb this sign into an in 
shift in Xp- 

Strictly speaking we have only given necessary conditions that some of the charge 
sectors in the tau-function satisfy (3.14), but extensive numerical investigation supports 
the conjecture that (3.32) and (3.36) provide a solution to (2.16). Moreover whilst we 
enforced fairly stringent conditions on a p to derive the necessary constraints on d p and C, 
many examples seem to confirm the idea that these values give a solution even when the 
conditions are not met. However, in those cases one generally might expect other solutions 
not of the form (3.32), (3.36). 

By analogy with the two-soliton case, we should proceed by discussing the position 
of the singularites in <p. Since the solutions corresponding to C± are parity conjugate, 
and we know from the negative energy that is singular somewhere on the real line 
— oo < x < oo, we know that one of the two solutions must be singular in the physical 
region x < 0. Unfortunately, finding the zeroes of the tau-function in general is very 
difficult. Numerical work suggests that it is more difficult to find regular solutions as more 
solitons are introduced. As alluded to above, solitons in the real coupling theory have 
negative energy, so that it seems that the lowest energy configuration consists of adding 
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in as many solitons as possible without making the solution singular for a; < 0. One 
possibility is that the best we can do is to place the singularity at the wall. In this case 
the analysis leading to the expression for the total energy (2.22) is still correct and yields 
a finite answer. Physically the infinities in the bulk and boundary energies exactly cancel. 

In parallel with the two-soliton case, one can ask which boundary condition does the 
solution given above correspond to; that is for what does the above solution satisfy 
the boundary condition. Again the answer seems to be a lot more complicated than for 
two-solitons, but nonetheless some of the features of that case do seem to persist. For the 
two-soliton case it was seen that one could write Tj(x = 0) in the form LWj where all the 
j-dependence is in W. From this we deduced that 



e^Vo = (3.37) 



>.*/2| _ \W 3 ^W j+1 \ 
and that 

Aj = ±sign(^_ 1 M/- +1 ) (3.38) 

where the plus/minus sign corresponds to the two solutions C±. We postulate that we can 
write the tau-function (3.30) evaluated at x = in the form LW? where we can expand 
W as 

Wj= E W,..^cosQ> p {f + ^-}) (3.39) 
By comparing this with the tau-function with d p given in (3.36), we find 

Ca lt r 2 ...a N = LI Sill ( —- AWrdr ~ <T s a s \)) (3.40) 

r<s n + 1 

and 

L = 4TT— ^ TT --. (3.41) 

We have checked that this ansatz for Tj(0) is true up to six solitons. By combining 
(3.40,3.39,3.38) we can deduce Aj for any given values of x P > and use this to deduce a 
'phase' diagram for different boundary conditions as a function of the Xp- One such ex- 
ample is given in Figure 2, where we have taken n = 8, iV = 2, ai = l,a 2 = 3. In this 
figure we have plotted the lines in the parameter space (xi,%2) along which r,-(0) = for 
j=0,..,8. 

In the region marked / the boundary conditions are found to be 

(A , A ± , .., A 8 ) = (-1, 1, 1, -1, -1, -1, -1, 1, 1). (3.42) 

The boundary conditions for the other regions can be deduced by noting that Wj changes 
sign as the line Tj = is crossed. So by (3.38), A,_i and Aj +1 change sign when this line 
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1 2 3 4 5 6 

Z 2 



Figure 2: Boundary conditions for N = 2, n = 8, a± — 1, a 2 = 3 

is crossed. One can check that the resulting boundary conditions are compatible with the 
Z n symmetry which is realised as 

- *-(^TT) (3 - 43) 

j - 3 + 1 (3.44) 

Suppose that we choose \v so ^at we lie on one of the lines Tj — 0. Then the corresponding 
solution for is singular, and in fact 

e ai -* = +oo, (3.45) 
e aj+K(f> = 0, (3.46) 
e a ' -1 ^ = 0. (3.47) 

Thus along this line the coefficients Aj_i and Aj +1 are undetermined as we may have 
expected. It follows that a particularly rich variety of boundary conditions are allowed at 
the points where many such lines intersect. 

If n is odd we can consider static soliton solutions containing a single constituent soliton 
of type a N+ i = (n + l)/2. Reality of the tau-function constrains the associated parameter 
£jv + i to be real, or equivalently Xp — 0, 7r. Following much the same argument as in the 
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previous section we initially assume that the a p = n+ 1 — a p , p = 1, .., N are chosen in such a 
way that only one term contributes to the sector of maximal charge Q ma x = Y, p =i o, p + a N+ i 
and two terms contribute to sectors of charge Q p = Q ma x — a p , p = 1, .., N. This results 
in similar expressions for C± and d p , p = 1, ..,N except that the sums in (3.32), (3.36) 
run from 1,..,N + 1. It only remains to determine the value of djv+i- Only one term 
contributes to the charge Qn+i sector since as there is no soliton conjugate to the middle 
soliton p — N + 1, we cannot add a soliton-conjugate soliton pair to the term involving 
solitons of type p — 1, .., N. Instead the resulting equation is 

N 

= (C± - ]T m ap f - 1 = (m aN+1 ± l) 2 - 1 (3.48) 
P =i 

Since fn aN+1 = 2sin(7r/2) = 2, we see that only the C_ solution is allowed. The value of 
djv+i is unconstrained by this or any other equation, so we are free to place the middle 
soliton anywhere! Under a parity transformation, the energy and the values of d p , p = 
1,..,N remain unchanged. This is in agreement with the idea that the energy of the 
solution on the whole line is formally given by the sum of the half-line energies of the 
solution and its parity conjugate; that is in this case 

N+l N 

-2(n + 1)(C_ + C_) = -2(n + 1)(2 £ m ap - 2) = -2(n + 1)(2 ^ m flp + m ajv+1 ) (3.49) 

p=i p=i 

as we should have expected from the masses of the constituent solitons. 

4 Scattering solutions for affl Toda theories 

In the previous section we showed how one can construct static background configurations. 
In this section we show how to solve the classical linearised perturbation equations around 
this background, and explicitly calculate the classical scattering matrix. Once more we 
shall rely on the tau-functions given by multi-soliton solutions. 

Let us assume that the field is infinitesimally perturbed about the vacuum configuration 

<j>(x,t) = <t>{x) vac + e(x,t). (4.1) 

Substituting this into the equations of motion and the boundary conditions (2.5), and 
keeping terms linear in e(x,t) yields 

r 

d^e(x, t) + J2 n i a i e ai -' t '—^ (a 1 ■ e(x, t)) = (4.2) 

i=0 
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and 

d x e(x,t) ^\Y j A i a i ^ i e ai ^^l 2 {a i ■ e(x,t)). (4.3) 

1 i=0 

Far away from the boundary x — 0, the field 4> vac {x, t) — >• for finite energy configurations, 
and in this limit the equation (4.2) reduces to (2.7) and the solution for e(x,t) tends to 

e(x, i) - p a e^(e^ + K a e"^) (4.4) 

where E 2 — p 2 = m 2 . This solution consists of the superposition of two oscillatory solu- 
tions corresponding to incoming and outgoing 'waves' or quantum mechanically 'particles'. 
On the other hand, the field <j>(x, t) corresponding to a single soliton of type a has the 
asymptotic form 

(f>(x,t) ~ Pa e*W (4.5) 

where $(x, t) = a(x — vt), and a 2 — a 2 v 2 = m 2 . Comparison of the two asymptotic 
behaviours suggests that we can obtain appropriate oscillatory solutions from the soliton 
tau-functions if we make the identification 

a = ±ip (4.6) 
av = %E. (4.7) 

From here it is clear how to proceed. We take the tau-function r vac corresponding to the 
static vacuum solution that we found in the previous section and add in two further time- 
dependent solitons with the identifications (4.7). If we label the incoming and outgoing 
solitons by indices J, and O respectively, then we take 

07 = —do = ip , ajVj = o vo = iE (4.8) 
dj = d = e , xi = -Xo = i> , a i = a Q = b (4.9) 

Note that these two solitons are not a conjugate pair but are both of species b. The relative 
phase between the two waves is given by 

K b = e- 2 ^ (4.10) 

This tau-function describes the scattering of a particle of species b on the background 
given by r vac . As described in section two, requiring that this tau-function satisfies one of 
the boundary conditions (2.5) up to 0(e) amounts to the equation 

n - t'I + 2Cr[ - [C 2 - 1)7-^=0 = 0. (4.11) 

Considering the cases where the vacuum solution has even/odd number of constituent 
solitons in turn, let us assume first that it has IN solitons. The linearity of equation (4.11 
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implies that we can divide Tj into charge sectors as in the previous equation, and use the 
equation (3.14) with the addition of a time- derivative term as in(4.11), i.e. 

f k - T' k ' + 2CT' k - (C 2 - l)T k \ x=0 = 0. (4.12) 

There are at least two highest charge terms in the tau-function with charge Q ma x = 
J2p=i cip + b of the form 

T Qmax = rU n ^'e'^^^M 1 

j=l l<j<k<N 



N N 



X 



e if e iEt+ipx -Q A ( aj b)^ + e -if e iEt-ipx -Q A (a 3 b)(_ p ^ j 

i=i i=i 



where A^ ajb \p) are calculated using the definition (3.7), (4.7). Substituting this expression 
into (4.12) we find that the scattering matrices on the two solutions corresponding to C± 
are given by 

A , = = nU^'iP) u 14) 

* 2, P ±miuU A< °''' ) i-py 

If the vacuum solution has an odd number of solitons, then we find that the scattering 
matrix is given essentially by the above expression for K_, but where the product runs from 
1, .., N + 1. Note that the reflection factor given by (4.14) only depends on the number and 
species of solitons in the background solution, not on the topological charge parameters \p- 
It follows that boundary conditions related by the Z n symmetry (3.44) will have identical 
scattering matrices. 

4.1 Classical reflection bootstrap equations 

In the introduction, it was pointed out that in two dimensions, integrability placed strong 
constraints on the S'-matrix of the theory, and that S and K satisfied various algebraic 
constraints. In the case at hand both S and K are diagonal and in the classical limit S 
tends to unity, and these two facts ensure that most of the algebraic constraints are satisfied 
automatically. However, one non-trivial check is the reflection bootstrap equation, which 
in the classical limit becomes 

K C (6 C ) = K a (6 c + i9 b ac )K\6 c - i6 a bc ) (4.15) 

where the fusion angles are given by 

, 7r(a + j>) ( , 

9 - ~ WTT) (416) 

6 = n/2-6 (4.17) 
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and a + b + c = n + 1. (For simplicity we ignore the case that a + b + c = 2(n + l). For more 
details about fusion angles see for instance [5] and about the classical reflection equation 
see [13,8].) 

In the equation (4.15), the argument of K c is the usual rapidity variable <fi which is 
related to the momentum p and energy E of the incoming particle of type a by the formulae 

p = m a sinh(0) (4.18) 
E = m a cosh(0). (4.19) 

Note that using the identification (4.7), we see that 9 P = <p + i| where 6 p appears in(3.5). 
If, as is usual, one defines a bracket notation 

sinhf^ 



t x ) = — q — (4.20) 

blI11 H2 2(n+l)> 



then one can write in this notation 

AW)(p) (^ + a 3 -b)(^-a J +b) 
A^ h \-p) (2±i + a, + b) (2±i - a, - b) 



(4.21) 



where the left hand side of this equation is one of the factors appearing in the scattering 
matrix K± in equation (4.14). A straightforward calculation shows that each of these 
factors individually satisfy the reflection bootstrap equation. To prove that K b satisfies 
(4.15) it remains only to show that the remaining factor 

2i P ~ < = I (4.22) 

2ip + ml (n + 1- a b )(a b ) 

satisfies the bootstrap equation, and this is indeed the case. 

5 Other simply-laced algebras 

In this section we shall make a few remarks about extending the results found in the 
previous two sections to affine Toda theories based on the D and E series of algebras. 
The major technical difficulty in extending the results to these algebras is that there is no 
explicit form for tau-functions which correspond to more than two-soliton solutions, and 
even those solutions that are known are considerably more complicated than those of the 
theories [17]. Because of this, we shall restrict our search to finding which single soliton 
solutions [2] can be found which satisfy the boundary conditions (2.5). The general form 
for the corresponding tau-functions is 

Tj = l + d^dpe" 1 ^ + .... + 5^{d p ) n ^e n ^ x (5.1) 
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where p labels the species of soliton as before. The coefficients 8f have been explicitly 
calculated for all affine algebras on a case by case basis, and they are found to be real if 
and only if the minimal representation X p of the corresponding finite Lie algebra is real. 
Restricting ourselves to such real tau-functions, we substitute the tau-functions (5.1) case 
by case into the equation (2.16) and solve for C and d p . The results are presented below. 
The labelling of nodes on the Dynkin diagram is given in Figure 3. 

5.1 d?> 

The only solution is associated with the triality invariant soliton of species p — 2. The 
corresponding tau-functions are given by 

r = Ti = r 3 = r 4 = 1 + d 2 e^ x , and r 2 = 1 - Ad 2 e^ x + {d 2 ) 2 e 2 ^ x (5.2) 

and we find the solutions C = y/E/2, y/E/2 ± 3^2/2, ^6/2 ± y/2/2 and correspondingly 
d 2 = —1, (5 ± 3\/3)/(— 5 ± 3\/3),2 ± \/3. Only the solutions corresponding to — 1, (5 — 
3^3)/(-5 + 3 v / 3),2- are singularity free in the region x < 0. 

For d 2 — — l,ro — > as x — > 0, so that the only non vanishing component of A2 • — > 
—00 in this limit. Nonetheless the solution has finite energy, and obeys the boundary 
condition Aq = A\ = A3 = A4 = — 1 and A 2 is unspecified since e a2 '^ vanishes at the 
wall. For d 2 = 2 — \/3 we find that r 2 vanishes at the wall so that \ 2 ■ (ft — > 00. We 
find that A , Ai, A 3 , A 4 are unconstrained and A 2 = 1. Finally the solution with d 2 = 
(5 — 3y / 3)/(— 5 + 3\/3) is regular at the wall, and obeys boundary conditions with Ai = — 1 
for all i. 

5.2 4 1} 

Again the only solution is associated with the species p — 2. The corresponding tau- 
functions are given by 

r = ti = t 4 = r 5 = 1 + d 2 e 2x , and r 2 = r 3 = 1 - 2d 2 e 2x + (d 2 ) 2 e 4x . (5.3) 

We find the solutions C — (1 ± ^2 ± \/2) . The corresponding values of d 2 are (2 ± 

y/2 ± v^)/(2 =f y/2±y/2). The solutions which are regular in the region x < are those 

with C — (1 — y^2 ± \/2) . These solutions satisfies the boundary condition with A = Ai = 
A 4 = A 5 = -1 and A 2 = A 3 = =f1. 
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n-1 n 



O 



o^d — a 

1 2 3 



o 



n-1 



O^D — O 



n-3 n-2 







2 



O — O-^O — O-^D 

1 3 4 5 6 



2 



o^D — o^D — c^y 

7 6 5 4 3 1 







o 



O^D — O^D — O^D — O^D 

8 7 6 5 4 3 1 



Figure 3: Root Labels for simply-laced algebras 
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5.3 dg\n>5 

The solitons corresponding to the tip nodes do not yield any solutions. The other solitons 
correspond to roots a p with n p = 2 and their tau-functions are given by the formulae 

r = -Ti = l + d p e mpX (5.4) 
cos( (2i ~ 1)7rp ) 

n = 1 + 2 . 2{n - l) , d p e m - x + (d p ) 2 e 2m * x ,2<t<n-2 (5.5) 

cos (2(A) 

r n _ 1 = r n = l + (-iyd p e m " x (5.6) 

The value of the coefficient involving cosine takes at least two different values for different 
i, so inserting the difference of two such distinct tau-functions into (2.16) immediately 
yields that 2C = m p ± y/2. Now substituting r and r 2 into the same equation we find that 
m p = and that d p = 2 ± y/3. Given that 

m p = 2v^sin( ^ ) (5.7) 

we find that we must have n — 3x + 1, p — 2x. With these values (5.6) reduces to 

r = -n = r n _! = r n = l + d p e^ x (5.8) 
r 3l = r 3l+1 = l + 2d p e VEx + (d p ) 2 e 2VEx (5.9) 
r M+2 = 1 - Ad p e^ x + (d p ) 2 e 2VEx . (5.10) 

Only the solution with d p = 2 — \/3 is non-singular in the region x < 0. The solution is 
compatible with the boundary conditions A 3i+2 — — 1, and the other Ai are unconstrained. 

5.4 <£> 

In this case only the solitons of species p = 2, 4 have real tau-functions. Solutions which 
satisfy the boundary conditions can be found in both cases. If p — 2, then m 2 = ^6 + 2^ 

and solutions can be found with C = m 2 /2,m 2 /2 ± ^6-2^/2 and d 2 = -1, 5 ± 2>/6. 
Only the solution with d 2 = 5 — 2^ is non-singular for x < 0, and satisfies boundary 
conditions (2.5) with A 4 = —1 and the other A { = 1. 

If p = 4, then m 2 = y^6 — 2-\/3 and solutions can be found with C = m 2 /2,m 2 /2 ± 



6 + 2\^3/2 and rf 4 = —1, 5 ± 2\/6. The solutions with d A = —1, 5 — are regular for 
x < and satisfy the boundary conditions A = A 1 = A 6 = At = 1 and the other A's 
unspecified, and — 1 for all i respectively. 



22 



5.5 4 1} 

Only solitons of species p = 3 and = \/6 can satisfy the boundary conditions with 
C = 1713/2,1113 ± y/2/2. The corresponding values of c?3 are —1,2 ± and only the 
solutions with rf 3 = —1,2 — \f?> are free of singularities for x < 0. The solutions with 
(I3 = —1 does not seem to correspond to any sensible boundary condition of the type (2.5) 
but for <i 3 = 2 — we find that A 1 = A 2 = A 6 = — 1 with the other A; L unconstrained. 

5.6 4 1} 

By using the two tau-functions corresponding to rii = 2 and also r , one can use similar 
arguments to those used for d£\ n > 5 to deduce that m p = yf(6i). Since there are no 
single solitons with such a 'mass' we find a contradiction so that there are no single-soliton 
solutions. 



5.7 Scattering 

In section four we showed how by adding in two further solitons with imaginary momentum, 
we could describe classical scattering solutions on the static 'vacuum' configurations. The 
tau-functions had to obey the equation 

ft - r'l + 2Cn^ - {C 2 n\ - nAr^o = 0. (5.11) 

As explained above, multi-soliton tau-functions for algebras other than are very com- 
plicated in general. However, the exception to this rule is r , which has exactly the same 
form as for the case, except that the interaction coefficient A^ pq \6) defined in (3.7) 
must be generalised. For a definition of the generalisation X^ pq \6) and a discussion of 
the properties it enjoys please see reference [20]. Restricting ourselves to single soliton 
backgrounds, and to linear order in the perturbation parameter e we find that 

T = 1 + d r e mrX + ee^e-^+iP* + ee ~^ e -iEt-i P x + d ^ e -iEt+i P x+m r x x {rb) ^ 

+ d r ee-^e- iEt - ipx+mrX X {rb \-p). (5.12) 
Substituting this into (5.11) and looking at the term linear in e one recovers 

(5.13) 

((C-m r ) 2 -l)((C-m r ) 2 + 2ip(C-mr)+m2-l)x( r6 )(p)-(C 2 -l)(C 2 +2ipC + m2-l)) 

6 Conclusions 

We have seen that integrable boundary conditions (2.5) and the equations of motion com- 
bine neatly to yield a simple equation satisfied by the tau-functions of affine Toda theory. 
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In the case of a large class of static solutions were found, and moreover it proved 
relatively straightforward to extract the scattering data on these backgrounds. As a consis- 
tency check we showed that the scattering matrices satisfy the classical reflection boostrap 
equation. 

Still, many open questions remain even in this simple case. We only provided a plau- 
sible argument that the conjectured form of tau-function did indeed satisfy the boundary 
condition, and we did not rule out the possibility of other solutions, amongst which the true 
vacuum solution may lie . The solutions consisting of 2N solitons contained N unspecified 
parameters Xp (those with 2N + 1 solitons were specified by A^ + 1 parameters). These can 
be thought of as moduli for zero-modes, since the energy of the solutions is independent 
of Xp- We can change the boundary conditions satisfied by the solutions by varying these 
parameters, in essentially the same way that one varies the topological charge of solitons 
in the imaginary coupling theory. The 'phase' diagram which specifies which boundary 
conditions can be obtained from a particular family of solutions by varying Xp seems in 
general very complicated and this makes it difficult to isolate the true vacuum solution. 
Nonetheless, the scattering data is found to be independent of the topological charge, so 
all the resulting boundary conditions will share the same scattering matrix K. 

A related difficulty is to find which solutions are singularity free in the region x < 0. 
The general rule seems to be that as we add in more solitons the energy of the solution 
decreases rather than increases because of the reality of the coupling, but also the solution 
develops more singularities. The trick to finding the true vacuum is therefore to add in as 
many solitons as one can without introducing singularities into the physical region. It may 
be that the best one can do is to have a singularity at the wall itself which is still physically 
acceptable. This may also provide a mechanism for removing the zero-modes of the vacuum 
solution, since it is conceivable that varying the parameters Xp would inevitably move one 
or more of the singularities in into the physical region. 

Although the formalism extends to other simply-laced algebras, the application to these 
algebras is hampered by the relatively complicated technology for tau-functions other than 
those of the series. Whilst some partial results were obtained in these cases, it is clear 
that more powerful methods are required. Clearly one way forward is to try and adopt the 
methods of [23] which provide exact (if not explicit) formulae for multi-soliton solutions 
which are equally neat for any algebra. In their notation the tau-functions can be expressed 
as 

T . = (A^e-^^ge-^-^-^A^e-^-^ (6.1) 
and the condition (2.16) can be written neatly as 

(AjKCn, + Ei)g(Crij - E^Aj) = (6.2) 
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where g is an element of the group associated with the affine algebra. It would be interesting 
to see if this yields a more complete solution to the problem. 

Finally, it is worth re-emphasising that analytically continued solitons seem to play 
an important role in real coupling Toda theory on the half-line and that this provides a 
respectable home for them. The model seems to be on a firmer footing than imaginary cou- 
pling theories with their manifest unitarity problems, and the singularities that inevitably 
occur for the solitons in the real-coupling theories can be placed behind the boundary out 
of harms way. 
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